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Математическое моделирование упругого скольжения в контакте 
колесо — поверхность дороги” 


М. В. Чувейко, С. В. Носачёв 


Рассматривается проблема упругого скольжения в контакте колесо — поверхность дороги. Разработана ма- 
тематическая модель взаимодействия поверхностей. Представленная модель предполагает, что деформаци- 
онными процессами поверхности дороги можно пренебречь. Колесо моделируется в виде совокупности сек- 
торов упругого обода, закреплённых на упругом диске. На динамику модели основное влияние оказывают: 
сила взаимодействия секторов обода с центром колеса, сила межсекторного взаимодействия и сила контакт- 
ного взаимодействия с поверхностью. Проведены экспериментальные исследования с использованием мето- 
дов численного моделирования. Выявлено, что существует зависимость величины срывной нагрузки от дру- 
гих свойств колеса. Результаты работы представляют интерес для решения задач управления движением 
колеса с позиции предотвращения срыва контакта. Применённый целостный подход позволит выявить 
эмерджентность, оказывающую влияние на результирующую траекторию системы. 

Ключевые слова: система колесо — поверхность дороги, деформация, сила трения, упругое скольжение. 


Введение. Проблеме динамики взаимодействия колеса с полотном посвящены многочисленные 
исследования. При этом анализ динамической системы машины основывается на изучении взаи- 
модействия этой системы с динамической связью, формируемой трибосопряжением в области 
контакта. В известных исследованиях динамическая характеристика фрикционного узла рассмат- 
ривается как заданная: её параметры не зависят от упругих свойств колеса [1—4]. В отличие от 
этих работ в предлагаемом исследовании доказывается, что свойства формируемой фрикционной 
связи не являются инвариантными к упругим характеристикам колеса. Тем самым уточняются и 
дополняются многие известные свойства системы. 

Данная работа посвящена уточнению закономерностей формирования упругих характери- 

стик обода колеса и их влияния на трибологические свойства пятна контакта, обусловливающие 
режим его функционирования (сцепление/скольжение). 
1. Построение математической модели. Рассмотрим колесо, представляющее собой обод, 
закреплённый на диске (рис. 1). Подобные колёса применяются на определённых видах спецтех- 
ники. Будем считать, что обод и диск являются упругими. Очевидно, что для математического 
моделирования подобной системы необходимо учесть динамику тангенциальных деформаций 
диска и обода колеса. 





Рис. 1. Схема колеса: С — сектор обода колеса; Дс — длина дуги сектора обода колеса; 
ДА/— длина достаточно малого сектора обода колеса 


* 
Работа по теме № 2964 проводилась в рамках выполнения госзадания. 
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Однако в случае малых смещений собственно тангенциальные деформации диска не пред- 
ставляют интереса, а имеют значение лишь с точки зрения их влияния на динамику деформаци- 
онных движений обода колеса. В связи с этим влияние диска колеса на отдельные элементы его 
обода будет описываться некоторой силой Ё., зависящей лишь от динамического состояния само- 
го обода (его смещения). В дальнейшем будем называть её силой взаимодействия с центром ко- 
леса. Для анализа динамики деформационных процессов математическая модель должна учиты- 
вать взаимодействие отдельных секторов обода колеса друг с другом. Силу, с которой один сек- 
тор обода колеса действует на соседний, будем называть силой межсекторного взаимодействия 
Е». Кроме того, математическая модель должна учитывать взаимодействие колеса с поверхностью 
дороги. Полагаем, что деформационными процессами поверхности дороги можно пренебречь, 
поэтому данная сила будет зависеть лишь от динамического состояния самого колеса. Будем на- 
зывать её силой контактного взаимодействия 2. 

Рассмотрим описанные силы более подробно. 

1.1. Сила взаимодействия с центром колеса. Рассмотрим малый сектор обода колеса С. В 
результате тангенциальных деформаций диска сектор смещается относительно своего естествен- 
ного положения на величину /[. Это приводит к возникновению «возвратной» силы Р., пропорцио- 
нальной длине дуги сектора АС и величине его смещения: 
Е“ =-— К ЁАС, 
где А. — коэффициент пропорциональности, в дальнейшем называемый удельной тангенциальной 
упругостью диска, кг/(с*'м). 

Определим величину /[. Для этого сектор С будем рассматривать как однородный цилиндр 
с объёмной плотностью О’, кг/м?. Пренебрежение кривизной сектора допустимо ввиду его малых 
размеров. Координаты поперечных сечений, привязанных к цилиндру, обозначим через х. В ре- 
зультате продольной деформации эти сечения приобретают новые координаты х =Ё (х) (рис. 2). 


То есть математически данный процесс может быть описан с помощью функции Ё которая пред- 
ставляет собой отображение множества исходных координат подвижных сечений {х,\ на множе- 


ство итоговых координат {х,!. Из физических соображений очевидно, что данная функция явля- 


ется гладкой и гомеоморфной, а следовательно имеет обратную функцию /Ё". 





Рис. 2. Схема смещения сечений обода колеса 
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Рассмотрим некоторое произвольное сечение цилиндра неподвижной плоскостью © с ко- 
ординатой хи. Определим величину натяжения / в данной плоскости. Под величиной натяжения в 
плоскости О будем понимать величину смещения подвижного сечения цилиндра (относительно 
естественного положения), отображающегося в данную плоскость, т. е. имеющего в результате 
деформации координату х,. Иначе говоря, величина натяжения [ может быть найдена по форму- 
ле: 

— 
ие (1) 
Из рис. 2 видны следующие соотношения: 
- - —1 
а о (х,). 


С их учётом формула (1) приобретает вид: 
И КО: (2) 


Выполним замену переменной интегрирования: 


ХЕ ЕЖЕ Ра = дк. 


Хх 


Тогда (2) будет: 
Х> 1 
[= | ——_дх. 
17) ‹3) 
Введём в рассмотрение величину х , определяемую выражением: 
х АЕ" (Хх 
и (4) 
дх>0 ЛХ Ах—>0 АХ 
С учётом ранее приведённых свойств функции Авводимая величина х существует в облас- 
ти определения Ги может быть найдена с помощью выражения: 
1 


Е 


Хх 


(5) 


Исходя из физического смысла т, в дальнейшем будем называть её коэффициентом уп- 
лотнения материала. Так как рассматриваются только продольные деформации, очевидно, что 
величина объёмной плотности цилиндра О после деформации может быть найдена по формуле: 





О(х)=О*(х). (6) 
С учётом (5) и (6) выражение (3) примет вид 
и [бек - 5 [9 ()гх. (7) 
Или, домножив и разделив на площадь сечения цилиндра Н: 
17 1 
[= х) Нах = —М.. 8 
= ] @(х) ЕЛ, (8) 
где Мь — масса материала, пройденного через неподвижное сечение © в результате деформации, 
КГ. 
Исходя из сказанного, величина возвратной силы может быть определена по формуле: 
КА! 
Е” = М,. (9) 
НО 


Однако данная формула применима только для неподвижного колеса. Для колеса, вра- 
щающегося с угловой частотой ©@, формула будет иметь вид 
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К. А/ _ 
К 
Г -^иб (м.-мМ,), (10) 
где М. — масса, проходящая через неподвижное сечение © в случае отсутствия деформаций. 


М. в момент времени & может быть найдена по формуле: 


о 


где А — радиус колеса, м; & — начальный момент времени. При этом полагается, что в момент 
времени & деформации колеса отсутствуют. 
Диссипативная реакция со стороны упругого диска может быть найдена по формуле: 
Е. ее ам (11) 
НО (Е 
где р, — коэффициент пропорциональности, в дальнейшем называемый удельным тангенциаль- 
ным коэффициентом диссипации диска кг/(с*м). 
Введём переменную $, в дальнейшем называемую потоком. Она представляет собой вели- 
чину массы, проходящей через некоторое неподвижное сечение © в единицу времени, т. е. опре- 
деляется выражением: 





ам 
Е (12) 
СЕ 
С учётом введённого обозначения совокупная сила Р., действующая на сектор со стороны 


упругого диска, может быть найдена по формуле: 








[9 = > 
Е =-^ к [5-Э%+1,(5-5)|, (13) 
НО | 
где $ — заданная величина потока: 
_ @М > 
$ =—® =о<КНО. 14 
о (в) (14) 


1.2. Сила межсекторного взаимодействия. Кроме возвратной силы на сектор С действуют 
силы, обусловленные неоднородным распределением деформации вдоль обода колеса. Речь идёт 
о силах взаимодействия с соседними секторами. Упрощённо данное взаимодействие можно про- 
демонстрировать как упруго-диссипативное взаимодействие абсолютно жёстких элементов массы, 
В дальнейшем называемых материальными точками. Рассмотрим некоторый момент времени &. 
Произведём разбиение обода колеса на равные секторы протяжённостью Д/ (рис. 1). Будем пола- 
гать, что Л/ имеет довольно малую величину, и неоднородностью деформаций внутри сектора 
можно пренебречь. Очевидно, что позиционная составляющая силы, действующей на гю матери- 
альную точку со стороны пружины, может быть найдена по формуле: 

Е (1) = (41-^/), (15) 


/ 
где А/, — протяженность пружины (го сектора) в недеформированном состоянии, м; А — ко- 


эффициент пропорциональности, в дальнейшем называемый удельной продольной упруго- 
стью, (кг*м)/с”. 
С учётом (4) данное выражение может быть преобразовано: 
А (16) 
А/ т, 


[ Г] 


Соответственно, диссипативная составляющая определяется выражением: 
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а1_, 1 ат, 
т а” 
где №» — коэффициент пропорциональности, в дальнейшем называемый удельным продольным 
коэффициентом диссипации, (кг*м)/С. 
Тогда сила, действующая на -ю материальную точку А, (1) со стороны пружины, может 
быть найдена по формуле: 


Е” (/)=-В (17) 


ен въ ат, 
а совокупная сила 
Е(1)=Е, (1 )-Е (+1 +2 (). (18) 
С учётом этого уравнение динамики материальной точки принимает следующий вид 
ть - та в, ен 


+1 











1+1 




















ле (19) 
--= ®|5-5%+1 ($-5)|, 
ПО 
где 7 — масса материальной точки (кг), которая может быть найдена по формуле: 
т = ОНА/. (20) 
Или, раскрыв скобки: 
[9 > в 
—- РКА -ВА - р ЖЕ к[$-5%+1 (5-5) |5 
СЕ" г. у @) НО; 
(21) 
ах Г. 1 ат. 
Н— =К, —А— —— к | $ - 5+1 ($-5 
во т ° [ т Но | -® | 
Устремим А/к нулю. Тогда уравнение (21) примет вид 
4“х д 1 [1 ат 
И =КА —-- А, — | — К [$ 5+1 ($-5 22 
“ав га В - я Г] | -® ) К 
С учётом (20) $ можно представить в виде 
ах 
$ =ОН —. 23 
и (23) 
Тогда уравнение динамики (22) примет вид 
а О Г Ш _ 
—=А —--И — || - К | 5-5 +В [$ -5$]|. 
а 0/1 ЕЕ не *1 ых ] 
Или, введя вспомогательное обозначение /: 
ь р, 2 - к ноАи +" ЧЕ] 
СЕ 9х ^0/\ =? @) НО а (24) 
ам _ 
0 
СЕ 


1.3. Сила контактного взаимодействия. Элементы обода колеса, находящиеся в зацеплении с 
поверхностью дороги, кроме уже описанных сил испытывают также силу, обусловленную трением 
движущихся поверхностей. Основываясь на исследованиях, представленных в [5, 6], величину 
данной силы определяем физическими характеристиками взаимодействующих поверхностей, а 
также скоростью их скольжения относительно друг друга. Качественно характер зависимости 
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удельной силы контактного взаимодействия (на единицу площади поверхности) Г, от скорости 
иллюстрирует рис. 3. 
Е, 
1 и 
09 | Г | 
0,8 | ы \ 
0,7 ! \ 
6 { \ 
0,5 | г | 
047} | 
ОЗ! | | 
0,21 | ый 
0,1! 


| | | | | | 


0 2 = 6 8 10 12 14 16 =103 М 
Рис. 3. Зависимость силы контактного взаимодействия от скорости скольжения 


Таким образом, сила, действующая на отдельный элемент обода колеса протяжённостью 

А/со стороны поверхности дороги, может быть найдена по формуле: 
Е, =, (у), (25) 
где и — скорость относительного скольжения м/с; /,(у) — функция, описывающая зависимость 


удельной силы взаимодействия от скорости, кг/с”. 
Или с учётом того, что поверхность дороги является неподвижной, а также используя 


формулы (23) и (6): 
Е. = А (26) 


5тР 


где №. =^, [== 
т 


С учётом сказанного уравнение динамики деформационных процессов обода колеса при- 
мет вид 


и | = О, (О, 


> | а 1 
№ = 
оо, 
ОЕ 


где У(/) — функция безразмерной величины, принимающая бинарное значение (0 или 1) в зави- 
симости от того, находится ли рассматриваемая точка / обода колеса в соприкосновении с по- 
верхностью дороги или нет: 


У(/) = 


1, / вобласти соприкосновения, 
0, —/вне области соприкосновения. 


1.4. Уравнение динамики коэффициента уплотнения материала. Получим уравнение ди- 
намики для т. Для этого рассмотрим некоторый сектор колеса протяжённостью Д/ ограниченный 
двумя неподвижными секущими плоскостями ©; и О›. Несложно видеть, что динамика массы /7 
рассматриваемого сектора определяется следующим уравнением: 
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от =5:-5, (28) 
где $, — поток массы через первое сечение колеса, кг/с; 5> — поток массы через второе сечение 
колеса, кг/с. 
С учётом (20) выражение (28) принимает вид 
О. 95. 
а НО м" 


Устремим протяжённость сектора Л/к нулю. Тогда уравнение (29) будет иметь следующий 


(29) 





Нм 5, о, 5 НА но = 


ВИД: 

Ат 1 05 

ЕЕ (30) 

СЕ НО 0! 
1.5. Общее уравнение динамики колеса. Уравнение динамики колеса как целого может быть 
легко получено из второго закона Ньютона. При этом внешними силами по отношению к колесу 
будут силы контактного взаимодействия с поверхностью дороги, а также некоторая сила сопро- 
тивления, действующая на колесо Ё.,. Будем рассматривать плоское движение колеса массы т 
вдоль горизонтальной оси. Тогда уравнение динамики координаты центра колеса х будет иметь 
ВИД: 


ах 
—_—_ = у, 
аЁ 
[4-Е р 
а 
ОЕ т | 


где 2 — протяжённость пятна контакта, м. 
Учитывая всё ранее сказанное, в частности уравнения (27), (30), (31), динамика движуще- 
гося колеса будет описываться системой уравнений: 


о 2 %)- - НБ АМ +В мА 
@ 0х 20/1) НО а | 
ам бб. 
[т 
Дт 1 05 
&« 09а’ (32) 
ах _ 
С _ р 
‚ |; а -—Е. 
Е т 


где 5 — поток, кг/с; 5 — заданная величина потока, кг/с; Н — площадь сечения цилиндра, м’; 
О — исходная объёмная плотность цилиндра, кг/м?; К, — удельная тангенциальная упругость 
диска, кг/(с”‘м); р, — удельная тангенциальная диссипация диска, кг/(с‘м); № — удельная 
продольная упругость, (кг*м)/с*; А› — удельная продольная диссипация, (кг*м)/с; М— величи- 
на натяжения, кг; х — коэффициент уплотнения материала; У(/ ) — функция зацепления; х— 


горизонтальная координата центра колеса, м; и — поступательная скорость центра колеса, 
м/с; 7 — масса колеса, кг; №, (и) — удельная сила контактного взаимодействия, кг/с^; В — 


протяжённость пятна контакта, м; Е, — сила внешнего воздействия, (кг*м)/с”. 
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1.6. Компьютерное моделирование. Компьютерное моделирование продемонстрировало не- 
тривиальный характер поведения представленной модели. Так, например, обнаружено, что ско- 
рость фактического движения центра колеса зависит от величины внешней нагрузки в режиме 
сцепления. Обнаружена также зависимость предельной до срыва величины внешней силы сопро- 
тивления (в дальнейшем называемой срывной нагрузкой) сцепления от скорости качения и т. д. 
Для демонстрации первого утверждения рассмотрим зависимость, представленную на рис. 4. 
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Рис. 4. Распределение тангенциальной скорости элементов обода колеса по его поверхности 


Данный график демонстрирует распределение тангенциальной скорости ИУ элементов 


обода колеса по его поверхности. Подобный характер распределения скоростей типичен для ка- 
чения деформируемого колеса в режиме полного зацепления с поверхностью дороги [7, 8]. Также 
на графике продемонстрированы две величины: фактическая скорость движения центра колеса 
И. и заданная скорость И. , то есть скорость, определённая классической формулой: 


У, =©оА, 
где д — угловая частота вращения вала колеса; К — радиус колеса. 
Несложно заметить, что имеет место соотношение И. >И. . Это приводит к возникновению 


некоторой ненулевой величины «скольжения», характеризуемой разностью скоростей: 
АИ=И, И. 
Однако в отличие от классического скольжения поверхности колеса и дороги в зоне их 
взаимодействия остаются неподвижными относительно друг друга, т. е. имеет место зацепление. 
Данный факт наглядно демонстрируется равенством И. =И., имеющим место в зоне контакта. 


Причина возникновения некоторой ненулевой величины ДУ — деформационные процессы, проте- 
кающие в колесе. Согласно результатам компьютерного моделирования динамика деформацион- 
ных процессов имеет довольно сложный характер, обусловленный большим числом внешних фак- 
торов и характеристик динамической системы колесо — поверхность дороги. Однако, несмотря на 
этот факт, зависимость ДУот величины внешней нагрузки имеет характер, близкий к линейному. 
Данное утверждение иллюстрирует рис. 5. 
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Рис. 5. Зависимость скорости центра колеса от силы внешнего воздействия 


На графике отображены функции зависимости фактической скорости движения колеса от 

величины препятствующей силы РА, для случаев, когда колесо обладает жёсткостями К\!, К, Кз. 
Очевидно, что характеристика менее жёсткого колеса обладает меньшей жёсткостью. Кроме того, 
неочевидным результатом моделирования оказался тот факт, что величина предельной до срыва 
силы зависима от жёсткости колеса при прочих равных условиях. Это свидетельствует о недоста- 
точности знания свойств и характеристик элементарного контакта трения для прогнозирования 
поведения динамической системы в целом. 
Заключение. В работе представлена оригинальная математическая модель системы колесо — 
поверхность дороги. Данная модель описана системой дифференциальных уравнений в частных 
производных. Она позволяет учитывать упруго-диссипативные свойства колеса, а также особен- 
ности взаимодействия контактирующих поверхностей в зоне пятна контакта. Математическая мо- 
дель системы позволила исследовать деформационные процессы, протекающие в колесе во время 
его качения, и обнаружить обусловленные этим процессом явления — такие, как упругое сколь- 
жение, возникновение устойчивого колебательного процесса в определённых режимах работы. 
Кроме того, небезынтересна обнаруженная зависимость величины срывной нагрузки от других 
свойств колеса. Факт наличия подобной зависимости позволяет говорить о невозможности сепа- 
рабельного рассмотрения динамических характеристик процесса трения контактирующих поверх- 
ностей и прочих свойств динамической системы для оценивания и прогнозирования характера её 
поведения. 

Иначе говоря, интерполирование трибологических характеристик элементарного контакта 
взаимодействующих поверхностей на трибологические характеристики всей системы является 
допустимым лишь в предположении абсолютной жёсткости колеса. 

Характер зависимости свойств трибологических подсистем и, в частности, величины срыв- 
ной нагрузки от свойств других подсистем, вовлечённых в динамический процесс, может быть 
предметом дополнительных исследований и представлять интерес для решения задач управления 
движением колеса с позиции предотвращения срыва контакта. Согласно результатам, приведён- 
ным в статье, подобный «целостный» подход к анализу и проектированию трибологических сис- 
тем позволит выявить эмерджентность, оказывающую влияние на поведение системы, а следова- 
тельно, и на её результирующую траекторию. 
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